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The sound speeds of solar models that include element diffusion agree with helioseismological 
measurements to a rms discrepancy of better than 0.2% throughout almost the entire sun. Models 
that do not include diffusion, or in which the interior of the sun is assumed to be significantly mixed, 
are effectively ruled out by helioseismology. Standard solar models predict the measured properties 
of the sun more accurately than is required for applications involving solar neutrinos. 



For almost three decades, a discrepancy has existed 
between solar model predictions of neutrino fluxes and 
the rates observed in terrestrial experiments. In recent 
years, the combined results from four solar neutrino ex- 
periments have sharpened the discrepancy in ways that 
are independent of details of the solar models Q. This 
development is of broad interest since a modest extension 
of standard electroweak theory, in which neutrinos have 
small masses and lepton flavor is not conserved, leads to 
results in excellent agreement with experiments ||]. 

Since the implications of a discrepancy with the stan- 
dard electroweak model are of great importance, the 
question persists: Can the solar neutrino problems be 
"solved" (or at least alleviated) by changing the solar 
model? This question has led to a series of generally 
unsuccessful ad hoc "Non-Standard" solar models || in 
which large changes in the physics of the sun are hy- 
pothesized in order to lower the calculated rate of the 
8 B neutrino flux. Over the past two decades, the most 
often hypothesized change is some form of mixing of the 
solar material that reduces the central temperature and 
therefore the important 8 B neutrino flux ||-[|. Previous 
arguments that extensive mixing does not occur are theo- 
retical, including the fact that the required energy is five 
orders of magnitude larger than the total present rota- 
tional energy (|]§,[l(| . Most recently, Cumming and Hax- 
ton |n[] proposed a flow of 3 He, characterized by three 
free parameters, designed to mix the sun in such a way as 
to minimize the discrepancy between solar neutrino ob- 
servations and predictions. By adjusting the parameters, 
these authors are able to reduce the calculated 7 Be flux 
more than the 8 B flux, a result not achieved in previous 
Non-Standard solar models. 

The diagnostic power of helioseismology [Q has been 
improved recently through the development by Tomczyk 
et al. lOl of an instrument that measures with the same 



equipment the low- and intermediate-degree mode fre- 
quencies. By providing a consistent set of frequencies for 
the lowest-degree modes, which penetrate to the greatest 
depth in the sun, these data constrain the properties of 
the solar core more tightly than earlier measurements. 

In this letter, we compare the solar sound speed c in- 
ferred from the first year of data jy| with sound speeds 
computed from standard solar models used to predict so- 
lar neutrino fluxes and find a rms agreement better than 
0.2% over essentially the entire sun, with no adjustment 
of parameters. Since the deep solar interior behaves es- 
sentially as a fully ionized perfect gas, c 2 cx T/fi where 
T is temperature and fi is mean molecular weight; thus 
even tiny fractional errors in the model values of T or fi 
would produce measurable discrepancies in the precisely 
determined helioseismological sound speed 
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This remarkable agreement between standard predictions 
and helioseismological observations rules out solar mod- 
els with temperature or mean molecular weight profiles 
that differ significantly from standard profiles. The he- 
lioseismological data essentially rule out solar models 
in which deep mixing has occured (cf. | fl5|| ) and argue 
against unmixed models in which the subtle effect of par- 
ticle diffusion-selective sinking of heavier species in the 
sun's gravitational field-is not included. 

Figure ^ compares the sound speeds computed from 
three different solar models with the values inferred 
fll2| , [l4j from the helioseismological measurements. The 
1995 standard model of Bahcall and Pinsonneault (BP) 
Jl?j| , which includes helium and heavy element diffusion, 
is represented by the dotted line; the corresponding BP 
model without diffusion is represented by the dashed line. 
The dark line represents the best solar model which in- 
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FIG. 1. Comparison of sound speeds predicted by different 
standard solar models with the sound speeds measured by 
helioseismology. There are no free parameters in the models; 
the microphysics is successively improved by first including 
diffusion and then by using a more comprehensive equation 
of state. The figure shows the fractional difference, 5c/ c, be- 
tween the predicted model sound speed and the measured 
[ p^|Jl^ | solar values as a function of radial position in the sun 

tis the solar radius). The dashed line refers to a model 
in which diffusion is neglected and the dotted line was 
computed from a model Jl(| in which helium and heavy ele- 
ment diffusion are included. The dark line represents a model 
which includes recent improvements in the OPAL equation of 
state and opacities |l7L|l8|. 



eludes recent improvements in the OPAL equation 

of state and opacities, as well as helium and heavy ele- 
ment diffusion. For the OPAL EOS model, the rms dis- 
crepancy between predicted and measured sound speeds 
is 0.1% (which may be due partly to systematic uncer- 
tainics in the data analysis). 

In the outer parts of the sun, in the convective re- 
gion between 0.7R Q to O.95i?0 (where the measurements 
end), the No Diffusion and the 1995 Diffusion model have 
discrepancies as large as 0.5% (see Figure |]). The model 
with the Livermore equation of state fl8|] , OPAL EOS, 
fits the observations remarkably well in this region. We 
conclude, in agreement with the work of other authors 
p9| , that the OPAL (Livermore National Laboratory) 
equation of state provides a significant improvement in 
the description of the outer regions of the sun. 

The agreement between standard models and solar ob- 
servations is independent of the finer details of the solar 
model. The standard model of Christensen-Dalsgaard el 
al. pc[ , which is derived from an independent computer 
code with different descriptions of the microphysics, pre- 
dicts solar sound speeds that agree everywhere with the 
measured speeds to better than 0.2%. 

Figure |l| shows that the discrepancies with the No Dif- 
fusion model are as large as 1%. The mean squared dis- 



crepancy for the No Diffusion model is 22 times larger 
than for the best model with diffusion, OPAL EOS. If 
one supposed optimistically that the No Diffusion model 
were correct, one would have to explain why the diffu- 
sion model fits the data so much better. On the basis 
of Figure [l| we conclude that otherwise standard solar 
models that do not include diffusion, such as the model 
of Turck-Chieze and Lopez plj , are inconsistent with he- 
lioseismological observations. This conclusion is consis- 
tent with earlier inferences based upon comparisons with 
less complete helioseismological data [^2|,^2 15 , including 
the fact that the present-day surface helium abundance 
in a standard solar model agrees with observations only 
if diffusion is included pq| . 

Equation [l] and Figure |l| imply that any changes 5T/T 
from the standard model values of temperature must be 
almost exactly canceled by changes Sfi/ fj, in mean molec- 
ular weight. In the standard model, T and /i vary, respec- 
tively, by a factor of 53 and 43% over the entire range for 
which c has been measured and by 1.9 and 39% over the 
energy producing region. It would be a remarkable coin- 
cidence if nature chose T and \i profiles that individually 
differ markedly from the standard model but have the 
same ratio T/fi. Thus we expect that the fractional dif- 
ferences between the solar and the model temperature, 
5T/T, or mean molecular weights, S/i/fj,, are of similar 
magnitude to <5c 2 /c 2 , i.e. (using the larger rms error, 
0.002, for the solar interior), 



\5T/T\, \S(i/fi\ < 0.004. 



(2) 



How significant for solar neutrino studies is the agree- 
ment between observation and prediction that is shown 
in Figure [If? The calculated neutrino fluxes depend upon 
the central temperature of the solar model approximately 
as a power of the temperature, Flux oc T n , where for 
standard models the exponent n varies from n ~ —1.1 for 
the p — p neutrinos to n ~ +24 for the 8 B neutrinos [g3| . 
Similar temperature scalings are found for non-standard 
solar models p4| . Thus, maximum temperature differ- 
ences of ~ 0.2% would produce changes in the different 
neutrino fluxes of several percent or less, much less than 
required [1| to ameliorate the solar neutrino problems. 

shows that the "mixed" model of Cummings 
(CH) Jll| (illustrated in their Figure 1) is 
grossly inconsistent with the observed helioseismological 
measurements. The vertical scale of Figure || had to be 
expanded by a factor of 2.5 relative to Figure [l] in order 
to display the large discrepancies with observations for 
the mixed model. The discrepancies for the CH mixed 
model (dashed line in Figure^) range from +8% to —5%. 
Since /x in a standard solar model decreases monotoni- 
cally outward from the solar interior, the mixed model- 
with a constant value of /i- predicts too large values for 
the sound speed in the inner mixed region and too small 
values in the outer mixed region. The asymmetric form 
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and Haxton 
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of the discrepancies for the CH model is due to the com- 
petition between the assumed constant rescaling of the 
temperature in the BP No Diffusion model and the as- 
sumed mixing of the solar core (constant value of /z). We 
also show in Figure |^ the relatively tiny discrepancies 
found for the new standard model, OPAL EOS. 
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FIG. 2. Non-standard solar models compared with helio- 
seismology. This figure is similar to Figure i] except that the 
vertical scale is expanded. The dashed curve represents the 
sound speeds computed for the mixed solar model of Cum- 
ming and Haxton |Hj with 3 He mixing. The dotted line 
represents the sound speed for a solar model computed with 
the rate of the 3 He(a,7) 7 Be reaction set equal to zero. For 
comparison, we also include the results for the new standard 
model labeled OPAL EOS in Figure |. 

More generally, helioseismology rules out all solar mod- 
els with large amounts of interior mixing, unless finely- 
tuned compensating changes in the temperature are 
made. The mean molecular weight in the standard so- 
lar model with diffusion varies monotonically from 0.86 
in the deep interior to 0.62 at the outer region of nu- 
clear fusion (R = 0.25R&) to 0.60 near the solar sur- 
face. Any mixing model will cause fj, to be constant and 
equal to the average value in the mixed region. At the 
very least, the region in which nuclear fusion occurs must 
be mixed in order to affect significantly the calculated 
neutrino fluxes Unless almost precisely canceling 

temperature changes are assumed, solar models in which 
the nuclear burning region is mixed (R < O.25i?0) will 
give maximum differences, 5c, between the mixed and 
the standard model predictions, and hence between the 
mixed model predictions and the observations, of order 



5c 
c 



(j,— < fi > 



7% to 10%, 



(3) 



which is inconsistent with Figure [l]. 

Are the helioseismological measurements sensitive to 
the rates of the nuclear fusion reactions? In order to 



answer this question in its most extreme form, we have 
computed a model in which the cross section factor, S34, 
for the 3 He(ct,7) 7 Be reaction is artificially set equal to 
zero. The neutrino fluxes computed from this unrealistic 
model have been used || to set a lower limit on the al- 
lowed rate of solar neutrinos in the gallium experiments 
if the solar luminosity is currently powered by nuclear 
fusion reactions. Figure |^ shows that although the max- 
imum discrepancies (~ 1%) for the S34 = model are 
much smaller than for mixed models, they are still large 
compared to the differences between the standard model 
and helioseismological measurements. The mean squared 
discrepancy for the S34 = model is 19 times larger than 
for the standard OPAL EOS model. We conclude that 
the S34 = model is not compatible with helioseismo- 



logical observations (see also Ref. 25 1). 

Some nuclear parameters are important for solar neu- 
trino experiments but have negligible effects on the com- 
puted solar model values of the sound speed. For exam- 
ple, we computed a standard solar model in which we 
artificially decreased by a factor of two the crucial cross 
section factor, Sir, f° r the rare 7 Be(p, 7) 8 B reaction. The 
sound speeds computed for this radically different value 
of S17 differ by less than 1 part in 10 4 from the standard 
model values. 

Finally, we comment on the effects of the recent im- 
provements in opacity and equation of state Q on 
the predicted solar neutrino fluxes. Table [j] gives the neu- 
trino fluxes computed for a series of three different stan- 
dard solar models, all of which include helium and heavy 
element diffusion. The model labeled BP95 is from Jl6| ; 
the models labeled New Opac and OPAL EOS include, 
respectively, the improved opacities discussed in |T^ ] and 
the improved opacities plus the new OPAL equation of 
state discussed in Ea| . 

TABLE I. Neutrino Fluxes for Solar Models with Diffusion. 
All fluxes, except for 8 B and 17 F, are given in units of 10 10 
per cm _2 s _1 at the earth's surface. The S B and 17 F fluxes 
are in units of 10 6 per cm _2 s _1 . 
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pep 


7 Be 


8 B 


13 N 


15q 


17 F 


BP95 


5.91 


0.014 


0.515 


6.62 


0.062 


0.055 


6.48 


New Opac 


5.91 


0.014 


0.516 


6.62 


0.062 


0.055 


6.48 


OPAL EOS 


5.91 


0.014 


0.514 


6.60 


0.062 


0.054 


6.45 



The neutrino fluxes computed with the improved opac- 
ity and equation of state differ from the previously pub- 
lished values [|l6| by amounts that are negligible in solar 
neutrino calculations. The predicted event rate, for all 
three models, is 



CI Rate = 9.5^ 



SNU 



for the chlorine experiment and 

GaRate = 1371? SNU 



(4) 



(5) 
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for the gallium experiments. The only noticeable change 
in the predicted event rates for the chlorine and the gal- 
lium experiment is a 2% larger event rate for chlorine, 
which is due to a small improvement pc| in the calcula- 
tion of the neutrino absorption cross sections for 8 B. 

We conclude that the recent improvements in opacity 
and equation of state do not significantly affect the cal- 
culated neutrino fluxes, although they do result in sound 
speeds near the solar surface that are closer to the mea- 
sured helioseismological values (see Figure |l|) . The calcu- 
lations of standard solar models lead to predicted sound 
speeds that agree closely with the measured helioseis- 
mological values. We cannot rule out with mathemati- 
cal rigor the possibility |27) of constructing nonstandard 
models, consistent with quantum mechanics and with 
other stellar evolution observations, that are tuned to 
give the same sound speeds as the standard solar models. 
However, Ockham's razor suggests a strong preference for 
standard solar models. 
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